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A bosonic field theory is derived for the tunable edge magnetism at graphene zigzag edges. The 
derivation starts from an effective fermionic theory for the interacting graphene edge states, derived 
previously from a two-dimensional interacting tight-binding model for graphene. The essential 
feature of this effective model, which gives rise to the weak edge magnetism, is the momentum- 
dependent non-local electron-electron interaction. It is shown that this momentum-dependence may 
be treated by an extension of the bosonization technique, and leads to interactions of the bosonic 
fields. These interactions are reminiscent of a 4> 4 field theory. Focussing onto the regime close to 
the quantum phase transition between the ferromagnetic and the paramagnetic Luttinger liquid, 
a semiclassical interpretation of the interacting bosonic theory is given. Furthermore, it is argued 
that the universal critical behavior at the quantum phase transition between the paramagnetic and 
the ferromagnetic Luttinger liquid is governed by a small number of terms in this theory, which are 
accessible by quantum Monte-Carlo methods. 



I. INTRODUCTION 

Graphene is a two-dimensional hexagonal lattice of car- 
bon atoms and shows many interesting features, rang- 
ing from Klein tunneling to anomalous quantum Hall 
behaviour X Also effects of electron-electron interactions 
in graphene have recently started to attract much inter- 
est. Because of the strong confinement of the electron 
wave functions to only one layer of carbon atoms, the 
interactions are very strong, as compared to other quasi 
two-dimensional systems. Yet, the large Fermi velocity 
and the vanishing density of states at the charge neutral- 
ity point efficiently suppress the manifestation of inter- 
action effects in bulk graphene. 

At graphene zigzag edges, however, states with a small 
kinetic energy, so called edge states, give rise to a strongly 
enhanced local density of states. These states are very 
susceptible to the electron-electron interactions: the ki- 
netic energy of edge states is small so that the physics 
of the electrons in these exponentially localized states is 
dominated by electron-electron interaction. The latter 
drives the edges of graphene to a ferromagnetically or- 
dered state, a phenomenon known as edge magnetism. 
This magnetic ordering has been predicted theoretically 
on the basis of many different methods for interacting 
electrons, such as mean-field approximations* 2 ^ ab-initio 
calculations*^ quantum Monte Carlo simulations^ and 
bosonizationj££ Recently, also experimental evidence for 
edge magnetism was reportedi^^ 

Usually, the edge states in graphene arc considered 
to have negligibly small kinetic energies as compared to 
the typical bulk states. Indeed, edge states are exact 
zero energy modes of the pure nearest-neighbor 7r-band 
hopping Hamiltonian3 In other words, the bandwidth 
of the edge states is zero in this simplest model. How- 
ever, it has recently been shown that the edge state band- 
width, and therewith the Fermi velocity vf of these one- 
dimensional states, can be tuned over a wide range, and 
in many waysi 8 ! 11 ' 12 A sufficiently large edge state band- 



width may reduce and even completely suppress the edge 
magnetism, so that the strength and also the type of the 
edge magnetism may be tuned experimentally^ 

In Ref. [H, it was shown that (tunable) edge mag- 
netism can be described on the basis of an effective 
one-dimensional edge state model which has been de- 
rived by projecting the Hubbard Hamiltonian of the two- 
dimensional honeycomb lattice onto the edge states. It 
was recently rigorously proven^ that the exact ground 
state of this model is ferromagnetic as long as the band- 
width is below a certain critical value. 

In an exact diagonalization analysis of this edge state 
models it was shown that the essential features of this 
model are (i) a momentum-dependent interaction vertex 
and (ii) the complete absence of umklapp scattering. Es- 
sentially, (ii) allows for the existence of ferromagnetism 
in one dimension, which is usually forbidden by the Lieb- 
Mattis theoremJ^ However, (ii) alone only allows a 1st 
order phase transition from a paramagnetic Luttinger liq- 
uid, appearing in the limit of large edge state bandwidths, 
to a maximally polarized state as for smaller bandwidths. 
This happens because, once the system is beyond the 
Stoner point (vp % U/tt, where Vf is the Fermi veloc- 
ity and U is the Hubbard interaction on the honeycomb 
lattice), the interaction energy gain associated with in- 
creasing the spin-polarization is always higher than the 
corresponding kinetic energy penalty, independently of 
the actual spin-polarization. Thus, the system with fea- 
ture (ii) alone flows to full polarization immediately as it 
hits the Stoner instability. The momentum-dependence 
of the interactions (i) provides a mechanism to stop this 
flow: essentially the interaction strength is inversely pro- 
portional to the spin-polarization, so that the interac- 
tion energy gain becomes smaller as the spin-polarization 
is increased. Thus, the flow to higher polarizations is 
stopped at a certain spin-polarization and the system is 
stabilized in an itinerant weak ferromagnetic stated 

On the basis of the identification of these two im- 
portant edge state features (i) and (ii), a generalized 
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edge state model was proposed and analyzed by exact 
diagonalization in Ref. Il2l In this paper, an interact- 
ing bosonic field theory is derived from this generalized 
fermionic edge state model, thus overcoming the restric- 
tion to small systems of the exact diagonalization analy- 
sis. For this the well-known bosonization technique^ is 
adapted to the momentum-dependent interactions. The 
proper treatment of the momentum-dependence of the 
interaction vertex is essential for the analysis of the weak 
edge magnetism. The resulting bosonic action contains 
interactions of the boson fields, giving rise to an appeal- 
ing interpretation as a generalized Landau functional of 
the form am 2 + 6m 4 + with m the local magnetization. 
The actual bosonic theory seems to be very complicated, 
but it will be argued that only few terms contribute to 
the critical properties at the transition between the fer- 
romagnetic and the paramagnetic Luttinger liquid. The 
critical theory, which is proposed for this transition, is a 
1+1 dimensional classical field theory with a real action. 
As such, this critical bosonic theory may be simulated 
with Monte-Carlo methods and does not suffer from the 
infamous sign problem. 

The paper is organized as follows. In Sec. [TTl the gener- 
alized edge state model and a simple fermionic mean-field 
analysis of it is discussed, in order to gain an overview of 
the phase diagram and the relevant mechanisms. Section 
IIIII contains the derivation of the bosonic theory, which 
is then analyzed in Sec. IIV1 A summary and a critical 
discussion of the results may be found in Sec. [V] 

II. THE FERMIONIC MODEL 

The derivation of the bosonic field theory starts from 
the generalized edge state model, introduced in Ref. [l2|. 
The umklapp scattering is not allowed for the edge states 
considered here, so that the Hamiltonian consists of three 
terms 

H = H + H{ s + # 1 bs , (1) 

describing the kinetic energy, the forward scattering and 
the backscattering, respectively. The linearized kinetic 
energy reads 

tt/6 

H =V F ( rk ) C lra C kra, (2) 

r=R,L fc=-7r/6 

<T="U 

where vf is the Fermi velocity, k is the momentum along 
the zigzag edge in units of the graphene lattice constant 
a ~ 2.4A, r = R, L stands for right- and left-moving 
electrons, respectively, and Ckro annihilates an r-moving 
spin-cr electron with momentum k. When used in an 
equation, r = ±1 for R, L, respectively. Note that k = 
is defined separately for left- and right-movers (see Fig. 
[T]). The domain of the edge states ([— f , f] for each r = 
R, L) is restricted to one third of the total Brillouin zone 
[— 7r, 7r] (see Ref. [l2t ). For convenience the Fermi level is 



set to zero energy, i.e. kp = for left- and right- movers. 
Note that, due to the absence of umklapp scattering, this 
choice can be made without loss of generality. 

The interaction Hamiltonian Hi = + Hi s is most 
conveniently defined in fc-space. The forward scattering 
Hamiltonian, corresponding to c/2 and 54 processes^ (see 
Fig. [T]) reads 

rrfs \ s \ 1 nr qt qt' or' w 

n l — ^ °k+q k k'-q k' x 

r,r' k,k' ,q 

'■ C i+qrt C krt c i'-qr'l C k'r'i , (3) 

where the /c-space summations are restricted such that 
all momentum arguments of the electron operators are 
in the interval [—§,§]■ L is the length of the edge in 
units of the lattice constant a and the total strength of 
the interaction is expressed by a Hubbard U . The colons 
indicate normal ordering of the fermion operators with 
respect to the non-interacting Fermi sea. The factors 

SI = y/l-r^k (4) 

with Ti G [0, — ] parametrize the momentum-dependence. 
For Ti = 0, the interaction vertex is momentum- 
independent. This limit corresponds to the forward 
scattering in a usual one-dimensional Hubbard chainJ^ 
ri = - corresponds to the interaction vertex of edge 
states. Although this work is finally targeted at the edge 
state limit ri = — , it is convenient for the bosonization 
to keep Ti as a free parameter. 




FIG. 1. (Color online) Left: The linearized dispersion, de- 
scribed by Ho- The dashed line indicates the original cosine 
dispersion of direct model derived in Ref. @. For Ti > 0, the 
electrons at high energy are weakly interacting while the low- 
energy sector is strongly interacting. The momenta of the 
left- and right-movers are centered around their respective 
Fermi momenta. Right: The allowed interaction processes <?i 
(backscattering), <?2,<?4 (forward scattering). Umklapp pro- 
cesses are forbidden for the edge states (see text). 

The backscattering Hamiltonian, describing g\ pro- 
cesses (see Fig. [J), is given by 

rrbs \ \ ^ \ 1 qt Q—r q—v qt w 

~ ~L 2-J k+q°k D fe'-g°fc' X 

r k,k',q 

C k+q,r,t Ck '- r 'f C k' ~q,-r.X° k ' ' r ^' (^) 

where the additional parameter Ab s has been introduced 
in order to be able to tune the backscattering strength. 
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Abs = 1 corresponds to the physical SU(2) invariancc, 
which is obtained by a direct derivation from the two- 
dimensional Hubbard model.— 

The factors enter H^ s in a similar manner as they 
appear in H[ s . In both cases, for each electron operator 
participating in the interaction, the factors rescale 
the interaction strength: high-energy electrons suppress 
the interaction and low-energy electrons increase it. This 
means that exactly those processes which usually invali- 
date the bosonization technique beyond a certain inter- 
action strength, namely the ones involving electrons far 
from the Fermi level, are suppressed^ Indeed, it turns 
out that a proper treatment of the S^. factors in the 
bosonization results in a well-controlled theory also be- 
yond the Stoner point, i.e. for strong interactions. 

In order to demonstrate the significance of T\, I pro- 
ceed with a simple fermionic mean-field (fMF) analysis of 
Ho + H^r^- The parameter Ak describes the splitting of 
the Fermi points kp ra = raAk of r-moving spin-cr elec- 
trons (u = ±1 for f- and 4_-spins, respectively, when used 
in equations). Obviously, Ak is proportional to the spin- 
polarization 771 = ^Ak £ [—1,1]. The variational energy 
as a function of the spin-polarization m, corresponding 
to a fMF treatment, is 



E fM F(m) 



U 2 LT?tt 2 



3() 



5184 



(0) 



For ttvf > U, EfMF(w) is minimal for zero polarization. 
This is the paramagnetic Luttinger liquid phase for weak 
interactions. For irvp < U, however, the ground state 
has a non-zero spin-polarization m^O. Fig. [2] shows m 
as a function of the Fermi velocity vf- 
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Compared to the fMF prediction for the position of the 
Stoner instability, exact diagonalization gives a transition 
shifted to higher U.— A previous analysis, based on a 
combination of fMF and naive bosonization^ failed to ex- 
plain this shift as this analysis was valid only sufficiently 
far from the Stoner point. One goal of the extended 
bosonization, presented in the following, is to overcome 
this limitation. 

Note that in Ref. [l3| the critical properties at the fer- 
romagnetic transition differ from the one discussed here. 
This due to different assumptions about the fillin g^ as 
the kinetic energy (vf in this work and A in Ref. Il3l ) 
is changed; in the present work, a constant number of 
electrons in the edge states is assumed, while Ref. [H 
assumes that the total filling approaches a trivial limit 
(completely filled or completely empty edge states) as 
the system becomes non-magnetic. 



III. BOSONIZATION OF THE EDGE STATE 
MODEL 

In its usual form^ bosonization maps an interacting 
one-dimensional fermionic model onto an essentially non- 
interacting bosonic models However, two prerequisites 
are essential for this mapping: (i) the single particle spec- 
trum must be linear, and (ii) the interaction vertex of a 
process gi (see Fig. [T]) must be a constant. For the edge 
states considered here, (i) is a good approximation, as 
has been shown in Ref. Il2t However, the momentum 
dependence of the interaction, which conflicts with (ii), 
is an essential feature of edge states, responsible for the 
stabilization of the weak ferromagnetism. Thus, (ii) is 
a priori not valid and the usual bosonization technique 
must be adapted in order to be applicable to the inter- 
actions of the form ([3]). 

Before discussing the momentum-dependent bosoniza- 
tion, it is instructive to revisit the usual bosonization and 
its breakdown at the Stoner instability. 



Naive bosonization 



FIG. 2. (Color online) Spin-polarization m from the fermionic 
mean-field theory. The different m curves correspond to dif- 
ferent momentum dependencies Ti = 0, — . 

The size of the spin-polarization is controlled by the 
m 4 term, which is proportional to T\. For Fi = 0, the 
magnetization instantly jumps from zero to its maximum 
value as the fMF Stoner point ttvf = U is crossed. In the 
m = 1 state all up-spin edge states are occupied and all 
down-spin edge states arc empty Such a scenario where 
the band edges determine the physics can evidently not 
be treated within bosonization. Fortunately, the Ti > 
in actual graphene zigzag edge states gives rise to a con- 
trolled intermediate regime of weak itinerant ferromag- 
netism (0 < | m | < 1), in which bosonization is applicable. 



If the interaction-induced electron-hole excitations are 
restricted to a small region around the Fermi surface in 
which the interaction vertex does not vary much, the lat- 
ter may be approximated by a constant, i.e., ~ 1. 
Reducing the functional form of an interaction vertex 
to constants is the central idea of g-ology, which is well 
justified for small interactions U. The resulting bosonic 
theory H s + H c + H± describes an ordinary spinful Lut- 
tinger liquid, with 



H u = ±- I dx 

Z7T 



u u K v {d x 6 v f + -^{d x 4> v ? 



(7) 



and u v K v = vp, u v jK y = vp(l ± U/irvp) for v — 
c, s, respectively. The bosonic field operators <j>„ and 
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6 V obey the commutation relations [0„(x'), d x 9 v i (x)] = 
iit5 vv i8(x — x') and their mapping to the fermion op- 
erators Ck ra is given in Appendix [Aj Their physical in- 
terpretation becomes obvious by reexpressing the spatial 
derivatives of the field operators in terms of fermionic 



field operators ^ r<T (x) 






ird x (j) c (x) = 


TO 


(8) 


ird x 9 c (x) = 


^r$ CT (x)<i/w(x) 

TO 


(9) 


nd x 4> s (x) = 


ra 


(10) 


ird x 9 s (x) = 


^m^ (T (x)V; r(T (x). 


(11) 



Obviously, the derivatives of the bosonic fields corre- 
spond to densities (</>) and currents (9) of charge (y = c) 
and spin (y = s). 

The backscattering Hamiltonian H^ s [see Eq. 
translates, for 5£ ~ 1, to a bosonic sine- Gordon 

term {2\ hs U/ {2iri]) 2 } f dx cos(V$«j)s) which is responsi- 
ble for restoring SU(2) invariance in the Luttinger-liquid 
regime^. The exact diagonalization analysis in Ref. [l2l 
suggests that, apart from restoring SU(2) symmetry, the 
backscattering is of minor importance in the magnetic 
phase. However, it cannot be excluded that H^ s changes 
the critical properties at the phase transition. Neverthe- 
less, the general idea of the bosonic theory can be under- 
stood on the basis of the forward scattering interaction 
only. Therefore, the breakdown of the naive bosonic the- 
ory [Eq. 0] is discussed on the basis of the forward 
scattering only. 

This naive bosonic theory breaks down for U > U cl -i t . = 
ttvf where the term in the Hamiltonian ([?])■ controlling 
spin density fluctuations (u s / K s )(j)' 2 , becomes negative. 
The reason for this breakdown is that the particle-hole 
excitations are not restricted to small energies for such 
strong U. In a sense, an interaction U > f7 C rit. drives the 
excitations to higher and higher energies, i.e. to regions 
far from the Fermi level (compare also the fMF analysis 
in Sec. |H]). As this interaction-driven flow to high en- 
ergies is not limited in one dimension, the particle-hole 
excitations hit the band edge. This regime in which the 
band edge becomes important cannot be described by 
bosonization. To see the breakdown on a purely formal 
level it is sufficient to note that for U > f/ cr it. the bosonic 
Hamiltonian is not bounded from below. 

For actual edge states, however, the interaction be- 
comes effectively weaker at higher energies because of 
the factors ££. Due to this suppression of the interac- 
tion, the flow to higher energies is stopped before it hits 
the band edge, so that this regime can be described by 
bosonization in principle. However, since the approxima- 
tion of constant SL is not valid, the usual bosonization 
technique must be extended in order to be able to ac- 
count for the ^-dependence of the interaction vertex. 



As described in Ref. H, it is possible to go beyond the 
critical interaction ?7 cr it. by a combination of fermionic 
mean-field theory and bosonization: the mean-field so- 
lution predicts a splitting between the up-spin and the 
down-spin Fermi level, which, due to the momentum- 
dependence of the interaction, effectively reduces the 
interaction strength for electrons at these new (spin- 
dependent) Fermi levels. Thus, if the bosonization is per- 
formed on the basis of the fermionic mean-field theory, 
it turns out to be stable also for U > ?7 crit . In a sense, 
the interaction is separated into strong parts which may 
be treated within mean-field theory, and weak parts that 
are accessible via the usual bosonization technique. 

This combined method, used in Ref. H, has severe lim- 
itations. Most importantly, it is only valid sufficiently far 
from the critical point. Also, the approximation of con- 
stant interaction vertices within the bosonization part 
of the method is not controlled. Thus, a more sophisti- 
cated approach is used in this work, treating the full mo- 
mentum dependence within the bosonization technique 
without resorting to mean-field theory. It is shown that 
the momentum-dependence of the fermionic interactions 
translate to interactions of the bosonic fields. 



B. Boson interactions from forward scattering 

Usually, interaction terms of the forward scattering 
form i-'l^'4'rt'4'i'i'4'r'i may be bosonized by means of the 



bosonization identity ipl a 4>r 



The Fourier 



transform of _ff( s to real-space, however, reads 

rr' 

s r (-k)^] [s r (k)i>rt] [s r ' (-k)^] [s r '(fc)M 



(12) 



with S r (±k) = Vl=F rT±k and k = —id x the momentum 
operator. Eq. (fT2")) may be expanded in powers of Ti so 
that a generalized bosonization identity for 



(13) 



for general integers n, m is needed. This is derived in 
Appendix Essentially the bosonic version of Eq. (fT3")) 
gives rise to a (n + ni + l)th power of first derivatives of 
the boson fields, so that Eq. (fT2")) leads to interactions in 
the bosonized theory. 

In principle, the bosonic momentum-dependent for- 
ward scattering has infinitely many terms. However, 
almost all terms involve at least second order spatial 
derivatives of the fields. These higher field derivatives 
correspond to spatial fluctuations of densities and cur- 
rents. Their size can be estimated for the largest al- 
lowed momenta; it turns out that these terms are small. 
Furthermore, one is foremost interested in the first or- 
der derivatives, since these are directly connected to the 
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global order parameters. Most importantly, the field <p s 
is related to the local spin-polarization m 

(#) = ^=m, (14) 

where \m\ = 1 for maximal polarization. In order to see 
the structure of the bosonic theory more clearly it is con- 



venient to partition H[ s = fff ,(0) + with H f ° M 

containing only terms involving the first derivatives of 
the boson fields, while all terms containing higher order 

derivatives are collected in H^ 1 '. 

The details of the momentum-dependent bosonization 
procedure can be found in Appendix Interestingly, it 
turns out that all orders higher than r 2 in H^ ' cancel 
exactly, so that 



H 



fs,(0) 



1 

UT 2 
8tt 







U 







I A 



V2n 



A'4 



29 



12 ±12 



29 



12 j/2 



29 



12 i/2 



+ 29: 



12 i/2 



29 



/2 i/2 



(15) 



Eq. (|15[) is one of the main results of this paper. It shows 
how the strong momentum dependence of the interaction 
vertex (ri 7^ 0) gives rise to interactions in the bosonic 
theory. 

The infinitely many terms in cannot be written 

down in a closed form. Fortunately, it turns out that this 
is not needed for the further analysis. For the leading 
order in Ti one finds 

= ±J dx^l K' 2 - 0f] + 0(T?). (16) 

The higher orders in Y\ also involve higher orders in k. 
As there was no sign of short wavelength instabilities in 
the exact diagonalization study of this models, one may 
safely restrict k <C T^ 1 and drop the 0(Tf) terms. 

On the first sight, the negative sign in front of the 
</>" 2 term appears to lead to a Hamiltonian which is not 
bounded from below. However, in practice, as k cannot 
be larger than T^ 1 and all k integrals are regularized by 
this UV cutoff, this term is always smaller than the (f>' 2 
term in H^°\ Also, it creates no mexican-hat minimum 

at finite k. Thus, the (f)" 2 term in H^ 1 ' may safely be 
ignored. 

In the critical theory H^ 1 ' is needed to control the 
short wavelength fluctuations in the correlation function 
(<p' s 4>' s ) a t the critical point. This is illustrated on the 
basis of a simplified classical model for the <j> s field. Con- 
sider the Hamiltonian density 

H = r0' s 2 + a0f + 0' s 4 , (17) 

where 4> s (x) is assumed to be a real classical field in this 
paragraph (and only in this paragraph). At the mean- 
field critical point of H, the prefactor t of the cf>' 2 term 
becomes zero. In the Landau mean-field theory, where 
the interaction 0' 4 is neglected, the term acj)" 2 is then 
the only non-zero term at the critical point, controlling 



the behavior of the correlation function 

(0' s (x)0' s (O)-0' s (O)C(O))^ ass 'cx J dk k ^~V . 

(18) 

For t = 0, the correlation function (|18p is proportional 
to —y/a. Note the similarity of (|T8| to the correlation 
function of a standard </> model. This is a consequence 
of the assumption of a classical field in %. Of course, 
the <f> 4 model in one dimension does not give rise to a 
thermodynamic phase transition. Treating (f> s correctly 
as a bosonic field conjugate to 9 S leads to a theory that is 
structurally different from a classical <fi A model; because 
of the anisotropy in time and space it is not possible to 
reduce the orders of the spatial derivatives by one, as in 
the classical model (fT7|) . 

The higher order terms in Eq. (|16[) involve higher order 
spatial derivatives, adding fc 6 , k 8 ,... terms to the denom- 
inator in (fT5|) . These do not change the infrared singular- 
ities and are thus, in analogy to the usual argumentation 
in the theory of critical phenomena,— expected to be ir- 
relevant for the critical behavior. 

Note that </>" 2 is not the only term which controls the 
field fluctuations at the critical point. The sine-Gordon 
term derived in the next subsection has a similar effect. 



C. Boson interactions from backscattering 

In this section, the backscattering processes are 
bosonized. For a constant interaction vertex, backscat- 
tering gives rise to a sine-Gordon Hamiltonian for the 
boson fields. Here, the momentum-dependence of the in- 
teraction vertex is taken into account, which leads to a 
generalized sine-Gordon Hamiltonian. The Fourier trans- 
form of the backscattering Hamiltonian [Eq. to real 
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space reads 



H* s = \ hs Uj2 



dx 



S r (-k)^ rt ] [s- r (k)i>- rt ] [S- r (-k)ft_ rl ] [S r (k)$ r 



(19) 



The bosonization of (|T9|) proceeds along similar lines as 
the bosonization of the forward scattering Hamiltonian: 
the S r factors are expanded in powers of Ti and each 
term is bosonized separately. Unfortunately, however, 
there is no cancellation of higher T\ orders as in the 
forward scattering case. In exchange, a closed formula 
for the bosonic Hamiltonian of the momentum-dependent 
backscattering process can be given 



ba 



u 



E/*» E 



Cm c n2 c n3 c n4 



ni ,ri2 ,n3 ,ri4 



[{irTAY'e- 21 ^] [(trT^Y^e 2 ^] [HrTA)"^- 2 ^] [(-irT \d x ) n * e 2 ^] (20) 



where c n is the expansion coefficient of the power series of 
y/1 — x = J2 n °nX n - The bosonic fields (j> r a are commut- 
ing combinations of the four bosonic fields introduced in 
Eqs. ® - jnj 



V8 



(r(f>c + ra4> s - C - a9 s ). 



(21) 



For more details see also Appendix [A"l 



From Eq. ((20)) all terms of the bosonic backscattering 
can be written down systematically. It is instructive to 
arrange the terms in powers of Ti and to reinstate the 
original bosonic fields CiS and (9 C S . Up to fourth order 
in. Ti, one finds 



H 



bs 



2U 
{2irqf 



dx 



i - Vzr^ 



Q/2 , 



cosV8(f> s + 0(Tf). 



(22) 



Note that the terms of order Tf vanish exactly in Eq. 
([22|) . The terms of order rf and higher contain higher 
spatial derivatives of the bosonic fields. At the critical 
point, they are not needed to control the fluctuations of 
the field as this is already done by lower order terms. 
Thus, according to the arguments given at the end of 
the last subsection, these are expected to be small and 
irrelevant for the critical behavior. 



A. Classical interpretation 

As already shown in Ref. [lU, the backscattering has 
only little impact on the phase diagram, so that one may 
start by dropping H^ s . In this subsection, the classi- 
cal vacuum of the energy functional Ec[<p s ,<pc,0 s ,dc] = 
Hq + H[ s shall be found, assuming all fields to be classical 
(i.e. [(/)„, V ] = 0). Only spatial derivatives of the fields 
enter the energy functional. Since the leading terms for 
the three fields 4> c ,6s, 9c have a positive coefficient, these 
fields are constant in the classical minimum of Eq and 
drop out of the further analysis. Thus, it is sufficient to 
consider 



IV. QUALITATIVE ANALYSIS OF THE 
BOSONIC FIELD THEORY 



In this section, an interpretation of the bosonic theory 
is presented. The goal of this section is to convey an 
impression of the significance of the individual terms for 
the ferromagnetic phase diagram. It should be empha- 
sized that the following handwaving arguments are by no 
means exact. The actual solution of the full interacting 
field theory requires sophisticated methods and is beyond 
the scope of this work. 



E, 



C[<Ps\ = 



1 

2^ 



v F 



12 



8tt ' 



url 



(23) 

Because of the positive sign in front of the </>" 2 term, 
the optimal classical field is a linear function of x, i.e. 
4> s (x) = ~^/^ mx i with m the magnetization (cf. the 

fcrmionic mean-field), so that Ec[4> s ] reduces to a Lan- 
dau function for m 



C\m] on 



irvp — U 
36" 



2^2 



-m + 



UTjir 
5184 



■m 



(24) 
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which exactly reproduces the fermionic mean-field theory 
of the generalized edge-state model. Of course, the ac- 
tual theory is more complicated because of the quantum 
nature of the fields. Thus, the features which are not cor- 
rectly predicted by the fermionic mean-field theory must 
be a consequence of the quantum nature of the bosonic 
theory. For instance, the classical treatment is not able to 
explain the dependence of the critical point on I\ which 
was found in the exact diagonalization treatment. 



B. Quantum fluctuations around the saddle point 
solution 

In this subsection, the quantum nature of the bosonic 
fields is taken into account on the simplest level: by a 
variational technique, the interacting Hamiltonian is ap- 
proximated by the 'closest' non- interacting Hamiltonian. 
This treatment is equivalent to a bosonic mean-field the- 
ory (bMF) and to a treatment of the quadratic fluctua- 
tion around the saddle point solution. Before this is done 
for the complicated Hamiltonian Hq + H{ s , a much sim- 
plified model of only one pair of bosonic fields 0, 9 with 
[4>(x'), d x 0(x)] = iir8(x — x') and a single </>' 4 interaction 
is discussed. Consider the Hamiltonian density 



H (0) - uK8' 2 



■ 12 



K 



(25) 



with uK = 1 and jl£ = 1 — U/nvp = r. Up to a global 
prefactor, Eq. (|25p contains the spin sector terms of or- 
der Y\ and one of the Y\ terms, i.e. a</>' 4 from H^'^. 
The 4>' 4 term has been chosen because r is allowed to 
change its sign. Then the 4>' 4 bounds the Hamiltonian 
from below. Following the prescription in Appendix [B] 
the bMF approximation of Eq. (f25j) is given by 



H 



(0) 
bMF 



q/2 



6a<0' 2 )]0' 2 , 



(26) 



so that one may introduce renormalized parameters u* 
and K* with 



u*K* 



1 



6a (<// 2 ) 



(27) 



At this point the mechanism which is responsible for 
the dependence of the critical point on Ti becomes obvi- 
ous: compared to the fermionic mean-field theory, which 
is critical at r = 0, the bosonic mean- field theory is crit- 
ical at t = — 6a (<// 2 ), with a oc Tf. In the language of 
the renormalization group, r — > r* corresponds to a mass 
rcnormalization, which is known to be non- universal, i.e., 
it usually depends on microscopic details. Here, these mi- 
croscopic details enter the renormalized r* via the cor- 
relation function (</>' 2 ) in which the fields have the same 
spatial and temporal arguments. On a formal level, a 
clear sign for quantities which depend on microscopic de- 
tails is the appearance of the high energy cutoff r\. 



Next, the correlation 



A'2\ 



ing from the bosonic theory of Luttingcr liquids^* one 



shall be calculated. Com- 

16 



would be tempted to calculate it directly from the bMF 
Hamiltonian (J2BJ. This leads to (0' 2 ) (O) = A'*/2r; 2 and 
it will turn out that this way of calculating the correla- 
tion function is not correct: K* = t* ' diverges at the 
critical point, and so does (<fi' 2 ^°\ Thus, with this naive 



average, Eq. (|27|) becomes r* 2 = tt* 3 + 3a and does 
not yield a proper solution with t* = 0, i.e., the system 
would never be able to reach its critical point. The rea- 
son that ((f)' 2 )^ is not the correct average is that at the 
critical point it diverges trivially in the sense that the 
Hamiltonian term governing the field <f> becomes zero for 
t* == and the Hamiltonian is independent of <f>. 

If the dominant term controlling the field <fi vanishes, 
previously subdominant terms become dominant and 
must be taken into account. The question is now, which 
out of the many terms in Eqs. (fTS"]). (fTB"]) . and (|22j) is 
the relevant one. Following the analogy to a Ginzburg- 
Landau theory, the quadratic second derivative term in 
Eq. (fl7J| is one candidate 



H 



(i) _ 



112 



(28) 



The average (0' 2 ), calculated with U [ °l lF + U {1) reads 



iJ2 



> = ;/ 



dfce-" |fc| 



(29) 



and is finite for t* 
r* = becomes 



= 0. The position of the critical point 



3a 



(30) 



The sine- Gordon Hamiltonian is also able to control 
the fluctuations of <j> at the critical point. However, 
Hamiltonians involving trigonometric dependencies on a 
field are difficult to treat. In this work I only give a 
very simplified account of the sine- Gordon Hamiltonian 
by sketching one possible scenario based on the assump- 
tion that the conventional RG treatment of sine-Gordon 
terms in spinful Luttingcr liquids 1 ^ extends to the present 
model. 

Consider a bosonic Hamilton density 



H = u*K*9' 2 



u 



9i± 



3(^0), (31) 



where cf> and 9 are the spin sector fields of a spinful 
Luttingcr liquid. g±± is the coupling constant of the 
backscattering processes. Note that in the presence of 
SU(2) symmetry there is a definite relation between K* 
and gi±. The interaction cos(v / 8</>) can be treated by a 
renormalization group analysis, within which K* — > 1 
and gi± — > 0. It is important to note that only for 
K* — > 1 the spin-spin correlation functions are indepen- 
dent of the direction of the spin-quantization axis chosen. 
In other words, the RG flow brings the system back to 
SU(2) symmetry. This calculation is perturbative in the 
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cos term and thus it is only valid for sufficiently small 
gi±. However, with the help of the SU(2) symmetry it 
may be argued that the validity of the RG idea can be 
extended to stronger interactions: only for K* — > 1 the 
theory describes an SU(2) invariant system and thus a 
system, in which this symmetry is inherent for all inter- 
action strengths, should be correctly described by K* = 1 
and <7i_l = 0, independently of the initial interaction pa- 
rameters. 

For the correlation function (0' 2 ) = K* /2ri 2 from 
the Hamiltonian u*K*6' 2 + u*cf> a / K* , this means with 
the sine-Gordon renormalizcd K* = 1 that (4>' 2 ) = 
l/2?7 2 , which gives a critical point 



Note, however, that Eq. (|32j) is by no means the result 
of a controlled calculation. It should be interpreted as a 
sketch of a possible scenario for the critical behavior of 
weak edge magnetism. 

Eqs. (|30]) and (|32|) display two different possible sce- 
narios for the critical behavior of edge magnetism. At 
this stage, however, it is not clear which scenario is the 
relevant one or if the fluctuations at the critical point are 
controlled by the combined action of the second deriva- 
tive and the sign-Gordon term. Moreover, it is not clear 
which term in the Hamiltonian controls the critical prop- 
erties at the transition. These questions are beyond the 
scope of this paper and are to be addressed in subsequent 
works. 



V. SUMMARY AND DISCUSSION 

An interacting bosonic field theory has been derived 
from the fermionic model of weak edge magnetism. In 
this derivation it is crucial to properly account for the 
momentum dependence of the effective interaction vertex 
function, which finally leads to the bosonic interactions. 
This makes the bosonization mapping used here different 
from the conventional bosonization (see, e.g., Ref. [ill ), 
where it is usually assumed that the interaction processes 
51,2,3,4 have an approximately momentum-independent 
strength. 

One of the most striking features of the bosonic theory 
derived in this work is the extended regime of validity. 
The usual theory of Luttingcr liquids, applied naively to 
the edge states, breaks down for interactions U that are 
strong enough so that U/irvp > I. But only above this 
bound, the physics becomes nontrivial by developing a 
spin polarization. Formally the breakdown of the naive 
theory happens because the bosonic Hamiltonian is not 
bounded from below in this case. For edge states, the 
momentum-dependence of the interactions leads to </^ 4 
terms with a positive prefactor. These terms ensure that 
the bosonic Hamiltonian is bounded from below and thus 
restore the validity of the theory for U/ttvf > 1- 



It is instructive to shed some light on the physical 
picture behind this mechanism. For a similar argumen- 
tation which is more targeted on an exact diagonaliza- 
tion analysis of edge magnetism, see Ref. [l2|. Electron- 
electron interactions soften the Fermi level at zero tem- 
perature by exciting particle hole pairs to higher ener- 
gies. For weak interactions, the kinetic energy penalty 
for these particle hole excitations is not overcompen- 
sated by the interaction energy gain. For sufficiently 
strong and momentum-independent interactions, how- 
ever, overcompensation sets in and the electrons are ex- 
cited to higher and higher energies by the interaction. 
Usually in one dimension, this process is not stopped 
until the the excitations hit the band edge. Thus, the 
band edge becomes important in this scenario, since it 
ultimately stops the flow of the system to high energy 
excitations. In bosonization, however, only effects ef- 
fects in which the band edges play no role can be de- 
scribed properly. This means that bosonization breaks 
down as the particle-hole excitations hit the band edges. 
For momentum-independent interaction this happ ens ex- 
actly at the Stoncr point (see also Ref. Il2i ). The 
strong momentum dependence of the interaction vertex 
of edge states, which makes the interaction effectively 
weaker at high energies (see Fig. [T]), stops the excita- 
tions before they hit the band edge. This is how the 
non-standard effective interaction vertex of edge states 
saves the bosonization from breaking down in the regime 
of weak edge magnetism. Thus, the special feature 
of strongly momentum-dependent interactions of edge 
states is crucial for the applicability of bosonization to 
onc-dimcnsional magnetism. 

The bosonic theory, consisting essentially of the three 
Hamiltonians (|T5|) , (fTo) , and (|22p , seems to be extremely 
complicated at first glance. However, it paves the way to 
a separation of important and unimportant terms. The 
universal critical behavior is expected to depend only on 
a small number of qualitatively important terms in the 
bosonic field theory, and a goal of future investigations 
should be to identify those terms and to extract the criti- 
cal properties of the ferromagnetic transition between the 
Luttinger liquid and the weak ferromagnetism. For in- 
stance, it was demonstrated in this work on a mean-field 
level that the terms T(f>' 2 +(j)'^ are not sufficient to describe 
the critical point. One rather needs to include a term in 
the Hamiltonian which controls the spatial fluctuations 
of the magnetization at the critical point. On the other 
hand, the term 9'^ in Eq. (|T5|) . for example, is expected 
not to affect the magnetic properties qualitatively, since 
a positive 0' 2 term already limits the fluctuations of 0' c . 

A thorough analysis of the singularities of this theory 
will give rise to the identification the qualitatively impor- 
tant terms. From the qualitative analysis of the bosonic 
field theory presented in this work, one may already state 
the preliminary expectation that the critical properties of 
the tunable edge magnetism is described a theory of the 
form 

Wait. = 0? + rfis + <t>'s + <t>? + 9 cos (33) 
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where the first term controls the quantum fluctuations 
(remember that 4> s and 9 S are conjugate fields), the sec- 
ond and third terms are the Landau function part of the 
theory, and the last two terms control the field fluctua- 
tions at the critical point. On the basis of this 1+1 di- 
mensional classical field theory, the effect of other terms 
on the critical behavior may be studied systematically, 
but it is expected that the essential physics at the tran- 
sition is captured by H C rit.- 

Note that the fields entering Eq. (j3"3")l are real and 
6 S enters only quadratically. Thus, S can be integrated 
out exactly and the critical theory can be formulated as 
a classical field theory with a real action. This enables 
the investigation by Monte-Carlo methods, since theories 
which can be formulated with real actions have no sign 
problem. 

It is furthermore important to note that the role of 
SU(2) invariance it is not clear at this stage. In the tra- 
ditional bosonization approach to Luttinger liquids^ the 
sine-Gordon Hamiltonian, describing the backscattering 
processes, is responsible for restoring SU(2) invariance in 
the correlation functions. A perturbative rcnormaliza- 
tion group treatment leads to a non-interacting Hamilto- 
nian with renormalized parameters (K* = 1) in the spin 
sector. The question in which way these arguments hold 
in the bosonic field theory of edge magnetism, derived 
here, is beyond the scope of this work. 
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Appendix A: Bosonization details 



fields ip ra . The goal is then to translate the typical terms 

^ 1 ^ 2 C 9 x 4 V'Jt(a ; i)Vrt(^2)^(^3)^(^) ! (A2) 

of which the expanded i/f s is composed, to a bosonic 
language. The four different spatial coordinates x\, 24 
are used in order to keep track of the correspondence 
between the differential operators and the fermion fields, 
and to be able to bosonize these terms by a point splitting 
method. In the end, the limit x\, X2, £3, X4 — > x must be 
taken. 

Since this work is concerned with Hubbard interac- 
tions only, the quartic Fermion terms always consist of 
two spin-up and two spin-down operators. This simpli- 
fies the analysis as the quadratic terms may be bosonized 
separately for each spin species. For the backscattering 
interaction [Eq. ((SJ)] the same idea as described above is 
used. However, the details of the procedure differ slightly 
compared to the forward scattering interaction. Note also 
that, as far as the bosonization is concerned, it is always 
assumed that the thermodynamic limit L — > 00 is per- 
formed in the end, so that only the leading terms in L^ 1 
are to be kept. 

1. The bosonic fields 

Following Ref. [l6l the fermionic fields may be ex- 
pressed in terms of exponentiated bosonic fields 



7F 



-*7™(aO p -*7r»(aO 



(A3) 



where r — R,L labels right- and left-moving Fermions, 
respectively, and a labels the spin. The Klein factors in 
Eq. (|A3[) have been dropped, since they appear only in 
the standard combinations for which it is well known that 
they are irrelevant^ The bosonic (but non-Hermitian) 
fields Jrcr(x) are defined as 



A generalized bosonization technique for momentum- 
dependent electron-electron interactions of the form ([3]) 
or ([5]) is derived. In terms of the real space fields 



L x / 2 ^ fc e' lkx Ckra, the forward scattering 



Hamiltonian reads 
H{ S = UJ2 J dx 

rr' 



(Al) 



where S r (k) = y/l — rT\k and k = —id x is the momen- 
tum operator, r = R, L labels right- and left-moving 
fermions, respectively. Here and henceforth, normal or- 
der with respect to the non-interacting ground state is 
implicitly assumed. _ff[ s is expanded in powers of the 
momentum operators, acting on the different fermionic 



where 5N r 



rSN ra Tr 
L 



EV^^V'' (A4) 

<7>0 ' q 



'fercr 



Ckra 



'kra 



Ckrc 



is the total 



density of r-movers with spin a, relative to the non- 
interacting ground state. The fc-space boson operators, 
in terms of the original fermionic operators Ckra , read 



Y4 Ckra . 

/ j k—rq,r,a Krcr 



(A5) 



Note that b qra is only defined for q > 0. On the ba- 
sis of the fermionic commutation rules for Cfc rcr , it can 



be showni^ that 



h 

Uqra, u n i r 'rr> 



At some 



Jqra, -q'r'a' 

points in the bosonization procedure, expressions need 
to be regularized. This is done in the standard wayi£ 
by introducing an exponential UV cutoff e - ''' 9 '/ 2 in the 



10 



/c-spacc sums. r\ has the meaning of a microscopic length 
scale, such as the lattice constant, thus a number of order 
one. In most of what follows, rj can be regarded as small 
compared to typical distances. However, if r\ needs to be 
specified explicitly, it should be chosen to be the inverse 
width of the reduced Brillouin zone, i.e. r\ ~ 6/77. 

In the final bosonized expressions, only the real parts 
of the four different 7 fields appear, i.e. 



1 



<t>ra{ x ) = 7;(lL( x ) + lra{x)). 



(A6) 



down-spin. The general form of the terms to be bosonized 
is 

F[n,m] = (id x ) n HcV) m : C(*)Vw(z') : (All) 

F[n, to] contains all factors corresponding to spin a in Eq. 
(|A2[) . The normal order : A := A — (A) of an operator 
A which is a quadratic form of fcrmion operators, with 
respect to the non-interacting ground state is needed to 
regularize the theory J£ The normal order arises directly 
from the mean-field treatment of the direct models From 
Eq. ([All) , one finds 



The spatial derivatives of rfa a may be interpreted as a 
local density of r-moving fermions with spin a, as is easily 
seen by substituting Eqs. (|A"4|) and (|A1)]) in Eq. (|A"6j) . 
The symmetric and antisymmetric combinations of left- 
and right-moving fermions give rise to the charge and 
current densities with spin a 

fa(x) = Y J {^{x)) 9 a {x) = - r{rfa a (x)). (A7) 



Furthermore, the symmetric and antisymmetric combi- 
nations of opposite spins form the charge and spin sector 
basis 



-±=(fa + fa) fa = -±=(fa - fa) (A8) 



V2 



V2 



(A9) 



Note that the field j ra does not commute with its com- 
plex conjugate jj, a . From Eq. (|A4[) . the commutation 
rules of the 7 fields may be calculated 



7r<r(2:),7rV'(y) 



Srr'Scra' log 



L 



277(77 — ir(x — x')) 

(A10) 



The subleading terms in L have been dropped. 



2. Forward scattering 

As explained above, the two factors for each spin 
species in Eq. (|A2[) may be treated separately The form 
of these two factors is the same for both, up-spin and 



il)l a (x)il) Ta (x') = 

cxpi(lL( x ') ~ 7^(^))expi(7 7 . g (x / ) - ~ir*{x)) . . 

2n(r] + ir(x' - x)) ' 1 ' 

where the two exponential factors in the middle have 
been commuted and thereby gave rise to the factor 
exp \y a {x)^l a (x')] = L/2ir(r) + ir(x' - x)) [see Eq. 
(|A10|) ]. After having applied the spatial derivatives 
{id x ) n (-id xl ) m to the right hand side of Eq. (jAT2l) . the 
expression is expanded into a Taylor series in 6x — x' — x. 
The terms of order Sx and higher vanish, the terms of or- 
der Sx° are the bosonic forms of F[n,m] and the terms 
of order 5x~ l , i > 1 are constant (and diverging) real 
numbers which are removed by the normal order. 

The terms of zeroth and first order in the derivatives, 
resulting from the procedure described above, are 



^[0,0] = 


r 

-' 

7T 






(A13) 


F[0,1] = 


-( 

7T 






(A14) 


F[1,0] = 


-( 

7T 




T 

+ id x —fa„ 


(A15) 



As expected, F[0, 0] is a local density term. The 
non-Hcrmitian terms, appearing in F [0, 1] and ^[1,0], 
are canceled in the final Hamiltonian. Note also that 
F[m,n — to] for different m but n fixed differ only by a 
total derivative (TD). This is enforced by the Hcrmitic- 
ity of the momentum operator —id x . For a Hamilto- 
nian of the form J dxF[n,m], these TD terms may be 
removed, as they only give rise to irrelevant boundary 
terms. However, since F[n,m] enters H[ s quadratically, 
the TD terms must be kept in principle. It turns out, 
however, that they do not give essential contributions to 
the physics. Therefore, we partition F[m,n — m] into 



F[m, n — m ] = F<°>[n] +d x F^[ 



m,n — to 



(A16) 



and discuss only the terms F^ for now. 
Each spin (er =^, I) contributes a factor 
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[S(ird x )4(x)] [S(-ird x ,)ti(x')] =J2J2(rFi) n c m c n - m F[m,n - m] 

r r n rn—O 

n 

= EE( rr i)" F(0) [ n ] E CmCn-m+d x ... (A17) 
r n m—0 



to the integrand of the Hamiltonian H[ s . c n = 2 n n\ 
is the prefactor of the nth order in the Taylor series 
of VT^x. The last expression indicates the TD terms 
which arc dropped from in the remainder of this section. 
For n > 1, the convolution of the c n gives zero, so that 

H i' (0) = ^ E / d * - ^ Tl<t> '^ ■ 

(A18) 

Transforming Eq. (|A18|) to the spin and charge sector 
basis, we arrive at Eq. (fFo)) , 

Appendix B: Derivation of the bosonic mean-field 
equations 

We use Feynman's variational principle in order to find 
the best approximation of an interacting theory by a non- 
interacting theory. More explicitly, the goal is to find a 
quadratic variational action S v for a given action S, so 
that the inequality 

pF = -hxZ = - In J V<Pe~ Sl 'e- {s ~ s >' ) 

= pF v -ln(e-( s - s ^ v <pF v + (S-S v ) v (Bl) 

is satisfied best. In other words, the goal is to find an ac- 
tion for which the right hand side of Eq. (IB1|) is minimal. 
Assuming that S = So + Si consists of a free part and 
an interacting part, we make the most general ansatz 

Sv =l%+ 2^BL E ^ ' M(fc ' w) ' ^ (B2) 



I 

where $>k,uj is a vector of spatial derivatives of fields in 
(k,u>) space. The matrix elements of M(fc,w) arc to be 
determined by 

S[/3F V + (S-S V ) V } = 0. (B3) 

As usual, this method gives rise to the following replace- 
ment rule for quartic interactions 

(B4) 

where the fields M arc components of the field vector $ 
and the contractions can be calculated by quadratic field 
averages 

Wr = (4f(z,T)<tr(x,T)) v . (B5) 

In the present case, the quartic interaction term is com- 
posed of <j)' , the spatial derivative of the bosonic field 
which describes the local spin density. Thus, all the six 
contractions are equal and one obtains 

atf* ^ 6a(</>' 2 )0 2 . (B6) 

Note that Eq. (|B4|) in its full generality can also be used 
for a mean-field treatment of the more complete model 
(fT5|) . It is, however, not obvious how the sine- Gordon 
term may be treated within this approach. All these 
issues are beyond the scope of this work and are to be 
investigated in future studies. 
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